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SUMMARY 

The natural -convection phenomenon is analyzed and it is found that 
the flow and heat transfer, in general, not only are functions of the 
Prandtl and Grashof numbers but also depend on a new dimensionless 
parameter. If this parameter is not negligibly small, the compression 
work and frictional heating may appreciably affect this mode of heat 
transfer , 

Consideration is given to the particular case of fully developed 
natural-convection flow of fluids with and without heat sources between 
two parallel long plane surfaces the temperatures of which are main- 
tained constant but not necessarily equal. These plates are oriented 
in the direction parallel to the generating body force. Solution of 
this problem yields detailed information on the velocity and tempera- 
ture distributions and heat transfer to be expected for such flows 
in tall narrow channels, on the effect of heat sources in the fluid, 
and on the effect of frictional heating on the process. It is found 
that the frictional heating and the heat sources increase the velocities 
and temperatures within the channel formed by the two surfaces. 
Increasing the ratio of the two wall-temperature differences (wall 
minus outside ambient) also leads to similar results. 


INTEODUCTION 

Flows which are generated entirely by the action of body forces 
(such as the gravitational force) on fluids with density variations 
due to heating are referred to as natural- or free -convection flows. 

It has previously been pointed out (see reference 1, for example) that 
natural-convection flows are of practical importance in aeronautics. 

The use of natural-convection flows in hollow passages in turbine rotor 
blades for cooling is one of the applications of this phenomenon in 
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practice. With the advent of the possibility of nuclear power, the 
natural-convection process becomes of even greater importance, because 
this mode of heat transfer appears in some of the many schemes for 
extracting the heat energy from an atomic pile. The use of liq.uid 
metals (in which heat may also be generated by heat sources) as the 
heat-transfer fluid for such applications is being considered because 
of their suitable behavior at the high-temperature levels that would 
be associated with atomic power. 

To date the theoretical investigations of natural-convection heat 
transfer have been restricted to such simple configurations as the 
single vertical flat plate and the horizontal cylinder. Further, 
the fluid considered in these investigations is usually air. The 
work done on more complex configurations, such as the natural- 
convection flow in channels or tubes, is for the most part experimental 
or semiempirical . An extensive experimental semiempirical study was 
performed by Elenbaas, who in reference 2 analyzed in an approximate 
manner the natural-convection heat transfer between two parallel plates 
heated to the same temperature and also made measurements for the case 
of the flow of air. In reference 3 semiempirical Nusselt numbers are 
compared for experiments of air flow in vertical tubes of different 
cross section, and in reference 4 (primarily a summary paper) 
functional equations for the Nusselt number are obtained by means of 
similarity considerations (that is, essentially by dimensional 
analysis). These equations are then rewritten by correlation with 
experimental or approximately computed results to yield semiempirical 
fonaulas for the Nusselt numbers. More recently, experimental 
investigations of natural-convection flows of liquid metals have been 
made; the results of such a study on horizontal cylinders are discussed 
in reference 5. None of this work on the more complex configurations 
yields detailed information on the velocity and temperature distribu- 
tions, and it does not apply for fluids containing internal heat 
sources. Further, the results predicted by these semiempirical 
formulas deviate in some cases from the existing experimental data. 

Therefore, in order to answer some of the many pending questions 
concerning the natural-convection flow of various fluids in enclosures 
and to obtain information on such flows of fluids containing heat 
sources, this phenomenon is analyzed herein. Particular consideration 
will be given to a simplified but representative case; namely, the 
natural-convection flow between two long parallel plates at consteint 
temperatures oriented parallel to the direction of the generating 
body force. This specific problem not only retains many of the 
physical characteristics associated with natural-convection heat 
transfer but also leads to a tractable mathematical problem. 
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ANALYSIS 

General Considerations 

In general, the differential equations governing the laminar steady 
flow of a viscous, compressihle, heat-conducting fluid which is subject 
to a body force are, in rectangular Cartesian tensor notation (see 
reference 6), 



li = p(T) 


(5) 


and 


k = k(T) 


( 6 ) 


(A complete list of the symbols used herein is presented in appendix A.) 
Equations (1), (2), and (3) e3q>ress, respectively, the conservation of 
mass, momentum, and energy; equation (4) represents a thermodynamic 
equation of state; and equations (5) and (6) represent the viscosity- 
temperature and thermal- conductivity - temperature variations. If p 
and k are assumed to be constant and if the coefficient of 
voluDietric expansion P is introduced in the body force tenn (see 
appendix B), equations (2) and (3) become 
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where Pd = P - P 3, e = T - T3, and the subscript s denotes a reference 
condition usually taken to be the hydrostatic condition. 


In a manner similar to that of reference 1 , eq.uations (2a) and (3a) 
(neglecting heat sources) can be written in dimensionless form as 



(3b) 


by letting Uj^ = Uj^/U, 0* = = Xj^/d and p = p/pU^. Here^ 

U should denote a unique velocity which characterizes the flow, and 
in the case of forced -convection flows this is taken to be the pre- 
scribed reference velocity Uc (as for example, the free-stream 
velocity) . From equations (2b) and (3b) it can then be seen that the 
solutions of the dynamic and thermodynamic problems for the forced- 
convection flow (at velocity Uoo), taking into account also the body 
force action, are given in terms of four parameters; namely, the 
Grashof number Gr, the Eeynolds number Be, the Prandtl number Pr, 
and the dimensionless temperature number 0 . The parameter 0 is 
defined 

^ , fa 

■ <=p(T„ - Tj) - 


and it is argued that if <<Cp(T^ - Tg) or^ eq,uivalently, if 

0^ < then the compression work and frictional dissipation terms 

(the last two terms, respectively, in ec[uation (3h)) can he neglected 
with respect to the conduction and convection. For the case of pure 
natural-convection flows, there exists no unique prescribed character- 
istic velocity Uoo, 30 that the parameter © as well as the Eeynolds 
number Ee becomes meaningless. Equations (l), (2a), and (3a) still 
hold, but now instead of using U^o as a reference velocity to obtain 
dimensionless equations, some group of the physical quantities (such 
as f^, and 3) which are directly connected with the natural- 


convection phenomenon must be utilized. Such a group was found in 
references 1, 2, and 8 to be 




( 7 ) 
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(For convenience^ is here taken to be the negative of the X-component 

of the body force per unit mass.) Hence, with this expression for U 
as the reference velocity, the dimensionless forms of the momentum and 
energy eq.uations for the case of natural-convection flow are 


and 


2^ 

3xJ “ Gr ^ Gr 


a f^l 



^PD 

^x'j^ 


(2c) 


1 he* 1 ^ pfxdfaui/aui ^uA 2 

y ^Xj Pr Gr ^Xj^Xj ” ^ °p ^ 5xj Cp j^^Xj \^Xj dx^ j ~ 3 ySxj j 


Pfyd 

Setting Gr = K, equation (3c) can be written as 


i 

y ^Xj Pr Gr dxjdxj 


_ 


K 




““J 


Su^ /Su 


du. 




(3c) 



(3d) 


The new dimensionless group ^fjd/c^ which appears here was first 

encountered in the analysis of reference 9. It does not appear explicitly 
in that reference^ however^ since for the particular case considered it 
was approximated hy Pgd/P hy means of the state eq.uatlon for a gas^ 
and this last expression was negllgihly small. This group was also 
obtained in reference 5 by a formal dimensional analysis; but its 
physical significance and function were not discussed. By comparison 
of equations (2c) and (3d) with (2b) and (3b) it can be seen that for 
the case of pure natural- convection flow the Grashof number is _ 
analogous to the Reynolds number for forced flows and the factor K 
is analog£US to the dimensionless temperature number 0. In fact; 

Gr and K are; respectively; simply Re and 6 based on U as 
given in equation (7). Therefore; for -pure natural- convection flows 
the influence of the compression work and frictional dissipation terms 
in the energy equation should be determined by the parameter K and 
not by 9 . Previously; the effects of the compression work and 
frictional heating had been neglected in the natural -convection 
phenomenon on the basis of qualitative arguments showing that 0 was 
always small; in order to do so a guess had to be made of the value__of 
the reference velocity. Since it has here been demonstrated that K; 
and not ©; is critical in determining the influence of the compression 
work and frictional heating on the flow and heat transfer in the purely 
natural-convection process; it must be determined whether in actual 
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practice K is always negligibly small. If K is always very small, 
then, of course, there would be no need to extend the previous work in 
this respect. However, calculations show that even under relatively 
mild conditions, moderately large values of K are possible. VThen it 
is further realized that the three prime physical factors in the 
natural- convection process, f^, p, and 9^, appear in the numerator of 

the expression for K and that in the more recent applications these 
could easily be many times those values usually associated with this 
mode of heat transfer (for example, could be as much as 10^ g in a 

centrifugal field, and p and 0^ also could be much larger, in atomic 

energy applications), it becomes clear that in many practical cases the 
compression work and frictional heating will influence the natural- 
convection flow and heat transfer. 

An interesting and different characteristic of the natural- 
convection phenomenon becomes evident if the compression work and 
frictional heating are taken into consideration. In this case, for 
example, the frictional heating is added to the physically imposed 
heat and should act as a heat source in the fluid and, hence, tend to 
increase the flow velocities. 


Specific Problem 

To solve the system of equations governing the natural-convection 
flow and heat transfer would at best be a formidable task because of 
the nonlinearity of the equations and because of the interrelation 
of the equations of motion with the energy equation. The consideration 
of the compression work and frictional heating terms in the energy 
equation makes the problem still more complicated. Therefore, in order 
to obtain equations which are tractable mathematically, it is necessary 
to make some simplifying assumptions. In this report, therefore, 
consideration is given to a simplified configuration which leads to’ 
less complicated equations but which, nevertheless, retains the 
essential physical behavior of the natural-convection process. In 
this way detailed velocity and temperature distributions can be 
computed, and the effects of heat sources and of frictional heat can 
be studied. 

Fully developed flow between long parallel plates with constant 
wall temperature . - The simplified configuration to be studied is the 
fully developed laminar natural-convection flow between two long 
parallel plane surfaces or plates which are oriented in the direction 
of the generating body force (see fig. l) and which are open at the 
ends to the ambient fluid. The surface temperatures are constant, but 
one surface may be at a different temperature from the other. For 
such a configuration it is assumed that the velocity, as in the more 
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familiar Poiseuille flow case, and temperature depend only on the trans- 
verse coordinate Y. The simultaneous assumption of these two conditions 
implies that there is always a net heat flow to the walls and also that 
the transverse velocity component V vanishes identically. Under these 
conditions equation (l) is identically satisfied, and the system of 
equations (2a) and (3a) becomes 


where now the more familiar notation = U, Ug = V, ^ 

X2 = Y is used. The product pp in the coefficient ^Pf^/ki in 

equation (8) can he written as W for fluids, as is discussed in 
appendix B. For gases, the density p can always be combined with 
the absolute viscosity p to form the kinematic viscosity v, and 
V should be evaluated at some convenient and representative reference 
point (as, for example, at the average of the wall temperatures). 

By equation (9), Pj, is seen to be a function of X alone. 

Since U and 0 have been assumed to be functions of Y alone, it 
is evident from equation (8) that dPj,/dX must be a constant. Hence 

the pressure gradient dP/dX inside the channel differs from the 
hydrostatic pressure gradient by at most a constant, since 


However, the pressure difference required to accelerate the fluid 
from the hydrostatic to the fully developed condition and the pressure 
difference to decelerate it back to the hydrostatic condition must be 
finite. Therefore, since the channel is assumed very long, the 
pressure gradient inside the channel becomes equal to the hydrostatic 
pressure gradient, and equation (8) may be written as 



( 9 ) 


( 8 ) 


and 



( 10 ) 


dP _ ^^s *^^1) _ ^ 

dX ■ dX dX ^ dX 


+ + constant 


dPD dPg 



(8a) 
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At this point several interesting ohservations can he made concerning 
equations (8a) and (lO). First of all, in equation (8a) as in the forced 
Poiseuille flow equation, the inertia terms (left side of equation (3c)) 
vanish, hut now the driving (buoyancy) term (second term in equation (8a)) 
is a function of the transverse coordinate; whereas, in forced 
Poiseuille flow the driving term (pressure gradient) is a function of 
the longitudinal coordinate. 

If the frictional heating term (second in equation (7)) is neglected, 
the energy equation is independent of the velocity distribution, hut the 
equation of motion, which yields the velocity, is dependent on the tempera- 
ture. This state of affairs is the opposite of that occurring in the 
forced-convection Poiseuille flow. Also, the convection and compression 
work terms in the energy equation (3a) now vanish. Although these last 
effects are eliminated, the solution of this simplified problem should 
not only yield practical results for the natural- convection flows in 
tall-narrow channels hut should also show the effect of heat sources 
in the fluid (last term in equation (10)) and the effect of frictional 
heating (second term in equation (lO)) on the natural-convection 
process . 

The boundary conditions associated with this problem are as 
follows; the fluid must adhere to the walls of the channel (the no-slip 
condition of viscous fluids) or, mathematically, 

U(0) = U(d) = 0 

and the temperature of the fluid at the plate must equal the plate 
temperature, or 


and 


0(0) = - T3 H 


0(d) = - T„ = mOy, = 

' ' S Wq 


Let 


U = (U/E)u = (k/fxPPd^)u 




K Pr 


u 


Y = yd 


e = (k|i/f/p2p2^4)^ = 


and 


2670 


2670 


9 


2V MCA TK 2863 


Equations (8a) and (lO) 'become 

u " + T =0 (11) 

X' ' + (u’)^ + oK = 0 (12) 

2 — 

where a = Qd /k0-^Q and K = Pr K = Pr Gr and the primes denote 

differentiation with respect to y. 

Hie boundary conditions are 

u(0) = u(l) = 0 
t(0) = K 

and 

t( 1) = mK 

where m = 6 /q . 

Wf Wq 

It is interesting that in this particular problem it is necessary 
to have a priori information on three temperatures (Tg, , and 

essentially T^^) in order to determine the actual temperature distri- 
bution^ and that these appear explicitly only as the ratio of two 
temperature differences in the analysis. Combination of equations (ll) 
and (l2j to eliminate the dependent variable x yields 

uiv _ (u')2 _ oK = 0 (13) 

with the boundary conditions 

u(0) = u(l) = 0 

u"(0) = -K (14) 

and 

u» *(1) = -mK 

For convenience, the heat-source distribution is taken to he 
uniform, so that Q and hence a are constants. The method of 
solution to he described is in no way dependent upon this restriction; 
that is, Q and a could he functions of the independent variable, 
and the same method could he applied in principle. In order to solve 
the given boundary value problem (equations (13) and (14)) a method 
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of successive approximations is employed. The eq.uation for this techniq.ue 
can he written 

- [u’(r_i)]^ - oK = 0 (15) 

where r = 0 , 1 , . . . denotes the order of the approximation and where 
u'_^ = 0 hy definition. It should he noted that r = 0 yields the same 

equation as would apply when the frictional heating is neglected. By 
elementary quadratures the zeroth-order (r = O) approximation Uq is 

4 Si 

>^0 = E Zl 7T (16) 

i=l 


where 

ai = ^ (4m + a + 8 ) 
ag = -1 



a^ = a 

and from equation ( 11 ) 

U7) 

The next higher approximation (considering now frictional heating) is 

k2 , 

1^1 = + 576(2520) 


where 

a^ = 35a^ + 320ma + 328a + 792m^ + 1392m + 840 

ag = 0 

83 = -2(51a^ + 444ma + 564a + 1008m^ + 2352m + 1680) 
a^ = 105 (4m + a + 8 )^ 
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^5 

= -1008 (4m 

+ a + 

8) 



^6 

II 

1 

(D 

6ma + 

6a 4- 

8m^ 4- 

8m - 64) 

^7 

= 24 (a^ + 4ma + 44a + 

72m - 

72) 

^8 

= 18 (3a^ + 

12ma - 

• 28a 

4-'12m^ 

' - 24m 4- 12) 


= -80a ^m 4 


) 



^10 

CD 

II 





and 







2 

■'l ■ ^0 - 576(3530) ^ 

Equations (16) and (17) yield the velocity and temperature distributions, 
neglecting the frictional heating; and equations (18) and (19) hold for 
small but significant values of K and, hence, show the effect of con- 
sidering the frictional heating to a first approximation. To obtain 
solutions for somewhat larger K, it would be necessary to continue 
with the iterative procedure described previously to obtain the 
higher-order approximations. Such a scheme becomes extremely tedious, 
and, furthermore, the convergence of the method can be established by 
comparison with a direct numerical solution of the complete boundary- 
value problem (equations (13) and (14)). Some discussion of these 
numerical results relative to the solutions given by equations (16) 
to (19) will be presented herein, and a more complete and detailed 
account of the numerical results together with a detailed account of 
the numerical solution procedure used is in progress. Comparison of 
these solutions obtained from equations (18) and (19), in which the 
frictional heating is considered to a first approximation, to the 
numerical solutions, in which the frictional heating is completely 
accounted for, should indicate the range of applicability and accuracy 
of equations (18) and (19). 

EESULTS AND DISCUSSION 
Velocity and Temperature Distributions 

The relations between the actual and dimensionless velocities and 
temperatures as determined from the various transformations in the 
analysis are 
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( 20 ) 


( 21 ) 


where U and 0 denote the actual^ and t and u denote the dimen- 
sionless q.uantities. For a given heat-transfer fluid and configuration 
(as specified through K and m) and for a given heat-source intensity 
(as specified through a) the dimensionless velocity and temperature 
distributions in equations (20) and (2l) can he computed from equa- 
tions (16) to (19). These computations will he accurate within the 
limits of the method of solution; that 1s, for small K. For 
larger K, computations of u and can hest he determined hy 
direct numerical solution of equations (13) and (l4). The range of 
applicability of the solutions given hy equations (16) to (19) will be 
discussed more fully subsequently. Eepresentative velocity and 
temperature distributions were calculated for values of K = 0.5^ 

3.0^ and 10.0. 

The particular values of ratio of wall-temperature differences m 
chosen for the computations correspond to the following interesting 
cases: (1) m = -1, in which the arithmetic average of the wall 

temperatures is equal to the temperature of the ambient fluid (that 
is, (T^^ + T^^)/2 = Tg); (2) m = 0, in which one wall is at the 

reservoir or ambient temperature; (3) m = 1, in which both walls are 
at the same temperature, and, hence, the effect of the wall-temperature 
difference is eliminated; and (4) m = 2, in which the wall temperatures 
are unequal but both are maintained at higher or lower temperatures 
than the ambient. 

Combinations of K euid m used in the computations and the figiure 
in which the results are plotted are given in the following table: 


U 


rp% 


and 


K 


e _ 


Wr 


K 

0.5 

0.5 

0.5 

0.5 

3.0 

3.0 

3.0 

3.0 

10.0 

10 . ol 

10.0 

10.0 

Eatio of wall -temperature 
differences, m 

-1 

0 

1 

2 

-1 

0 

1 

2 

-1 

0 

1 

2 

Velocity-profile figure 
number 

2(a) 

2(b) 

2(c) 

2(d) 

3(a) 

3(b) 

3(c) 

3(d) 

4(a) 

4(b) 

4(c) 

4(d) 

Temperature -prof i le 
figure number 

5(a) 

5(b) 

5(c) 

5(d) 

6(a) 

6(b) 

6(c) 

6(d) 

7(a) 

7(b) 

7(c) 

7(d) 
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For each K and m comhination profiles were calculated for a = 0^ 
10, and 100 with frictional heating neglected (given hy equations (16) 
and (17) and denoted by uq and Tq on the figures), with frictional 

heating included to a first approximation (given by equations (18) 
and (19) and denoted by u^ and '^i), and, in several specific cases, 

with frictional heating completely accounted for (given by the 
numerical solution of equations (13) and (14) and denoted by u 
and T ) . 

The computed results presented herein pertain specifically to a 
configuration wherein the body force is acting in the negative 
X-direction and at least one wall temperature is always greater than 
the ambient (that is, T^>Tg). 

Effect of different wall-temperature configurations (m varying) 
and heat sources (a varying) . - Examination of typical cases in 
figures 2 to 7 shows that, as expected, an increase in the heat-source 
parameter a or an increase in the wall-temperature parameter m 
results in larger velocities and higher temperatures. The velocity 
profiles change in such a way that the net mass throu^-flow, as 
represented essentially by the area under the u-curves, increases 
with m and a from zero at m = -1 and a = 0. The velocity 
profiles become more symmetric with increasing a, and for a = 0, 
negligible frictional heating, and m = 1, the profiles are similar 
to the Poiseuille profiles. For any given set of conditions, a 
decrease in the net mass through-flow or even no net mass through-flow 
can, of course, be obtained by adjusting the wall -temperature ratios 
so that m takes on larger negative values. It can also be seen 
from the temperature distributions that, as previously predicted, in 
all cases considered, either heat is being transferred to both walls 
or heat is flowing from one wall out through the other. For the case 
of no heat sources in the fluids and neglecting frictional heating, 
the temperature distributions, as are to be expected, are just the 
conduction profiles. 

Effect of frictional heating . - The effect of frictional heating 
can be seen in figures 2 to 7 by comparing the curves computed by 
neglecting frictional heating (denoted by subscript zero) with those 
computed by including frictional heating as a first approximation 
(denoted by subscript unity) for a given set of conditions. Numerical 
solutions of equations (13) and (14) obtained for several specific 
sets of conditions are also included (with no subscripts) for com- 
parison with the approximate solutions. For K = 0.5 and ail a 
and m considered, the frictional heating effect is small; that is, 
the U]^ and and u and x curves are not appreciably differ- 

ent from the Uq and Xq curves, respectively. However, even for 


14 


MCA TN 2863 


values of K relatively near unity, as is here represented by the case 
K = 3, there are conditions of m and a in which the frictional 
heating begins to alter the results appreciably. This situation could, 
depending on m and a, also occur for even lower values of K. In 
fact there are many combinations of K, m, and a for which the 
frictional heating becomes important. For example, the deviation 
between the solutions neglecting frictional heating (denoted by subscript 
zero) and those including the frictional heating to a first approxi- 
mation (denoted by subscript unity) for the case a = 0 first become 
marked for K = 10 and m = 1 (see figs. 4(c) and 7(c)). This deviation 
is more pronounced the higher the value of K for given _m and a. As 
was previously stated, moderately large values of K (or K) can be 
obtained under relatively mild conditions. For example, for the natural- 
convection flow of air under the influence of gravity at room tempera- 
ture, with 9 ^ equal to 1000° R and the Grashof number equal to 10^, 

K is approximately 12.1 (K = 16.8); for the flow of water in a gravi- 
tational field with of 150° R and a Grashof number of 10®, K is 

approximately 26.1 (K = 3.7). From the specific computations made 
herein, it can be seen that for a = 10, the deviation mentioned previ- 
ously first becomes apparent for K = 3 and m = 2 (see figs. 3(d) and 
6(d)); and for a = 100, a large difference exists for K = 3 and m = -1 
(figs. 3(a) and 6(a)). Hence, on the basis of the computations presented 
herein, it can roughly be stated that the solutions as given by equa- 
tions (l6) and (17) will be sufficiently accurate up to the limits 
previously stated. Beyond these limits the solutions including frictional 
heating effects to a first order (equations (18) and (19)) should be used. 
It should be kept in mind, however, that more detailed computations are 
necessary to define these limits more precisely. 

Numerical solutions which completely include the frictional heating 
effects and which are denoted by u and T without subscripts are also 
presented in the figures only for the cases of (K, m, a) of (3, -1, O) , 

(3, 1, 0), (10, -1, 0), (10, 1, 0), (10, 2, 0), (3, 1, 10), (3, 1, 100 ) , 
and (lO, 2, lO) in order to show the relative accuracy of the uj and 
solutions. For most of these cases, these numerical solutions 
coincide with the Uj^ and solutions; but in the higher internal- 

heat-generation cases (as given by high K, m, and a), these two sets 
of solutions differ (see figs. 3(c), 6(c), 4(d), and 7(d)). A signi- 
ficant difference indicates that the u^ and solutions, as given 

by equations (18) and (l9), respectively, are not sufficiently accurate. 
For the range of values of K, m, and a computed here, it appears 
that equations (18) and (19) could well be used for a = 0 for all 
K and m in this range, for a = 10 up to values of K = 10 and m = 2 
(figs. 4(d) and 7(d)), and for a = 100 up to K = 3 and m = 1 
(figs. 3(c) and 6(c)). Beyond these limits, equations (18) and (19) 
lose their accuracy, and the complete numerical solutions of equations (13) 
and (14) should be used. 
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The consideration of frictional heating not only changes the shape 
of the temperature profiles hut also leads to flovs with increased 
velocities^ as was expected. This change of temperature profile 
appreciably alters the temperature gradients at the walls and hence 
the heat transfer. 


In addition to these trends, the consideration of frictional heating 
leads to other interesting results. These results, which were obtained 
in computing the numerical solutions which completely take into account 
the frictional effects, are that two solutions exist for a given set 
of conditions and that there exists a critical set of conditions beyond 
which no solutions exist. Physically^ the presence of two solutions 
predicts the existence of two heat-transfer and flow states for a given 
set of conditions, which appears to be a situation analogous to that 
in a Laval nozzle. The existence of critical conditions appears to be 
similar to the thermal- choking phenomenon. One of the two flow states 
mentioned (the one encountered first physically) corresponds, as is 
seen in the figures, to the u^^ and solutions (or, more accurately, 

the u and t solutions); but the other, computed here only for the 

case of K = 10, m = 2, and a = 0, denoted by the curves labeled 
( 2 ) 

and T on figures 4(d) and 7(d), represents velocities and tempera- 
tures many times larger than the others at the same conditions. These 
results are being further investigated. 

Actual values of the velocities and temperatures can be computed 
from equations (20) and (2l) and figures 2 to 7. Since so many factors 
appear in those expressions, it is clear that there is considerable 
freedom in choosing the fluid and the physical conditions to obtain 
almost any given flow and heat transfer. For example, if Gr = 2X10®, 
fj = 6g, and 6^^ = 100° E for air at standard room conditions 

(Tg = 500° B) , a maximum velocity of 382 feet per second and a maximum 

temperature of almost 700° E would be obtained with both walls at the 
same temperature (that is, m = l). These conditions could very easily 
be obtained in an extremely mild centrifugal force field. 


For such computations to have quantitative significance, it should 
be kept in mind that this analysis pertains only to laminar flows. It 
is not possible to state under what conditions the transition from 
laminar to turbulent flow will occur, because no general stability 
theory exists for natural-convection flows. The only available natural- 
convection stability information is the experimental result indicated 
in reference 10 that transition on a single vertical plate occurs at a 
Grashof number of approximately 10®. (All values of Gr used in the 
numerical examples herein were less than 10®.) There is, in fact, 
relatively little general stability information for forced -convection 
channel flow with heat transfer, so that it is not even possible to 
obtain a rough transition criterion by relating the present problem to 
an equivalent forced-convection problem. 
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Although computations are presented herein only for a particular 
configuration, the analysis itself and the approximate solutions given 
hy eq.uations (16) to (19) are in no way limited to a single configura- 
tion. A change only in the body force direction (sign) merely alters 
the flow direction; that is, the sign of the velocity changes, or, in 
mathematical terminology, the velocity is an odd function of the body 
force. If there are no heat sources in the fluid and the effects of 
frictional heating are negligible, a change in the surface thermal 
condition (say, from T^>Tg to ‘I^<Tg) would also result in a 

change in flow direction. However, when either heat sources are 
present or frictional heating is not negligible, or both, the anti- 
symmetry with respect to 6 ^ is disrupted. For the case T^> Tg, the 

internal heat (due to heat sources or friction) increases the flow in 
a given direction as is shown in figures 2 to 4; but for '^Ts, this 
heat tends to retard the flow in a given direction and, if large enough, 
can change the flow direction. These effects can be seen in figure 8, 
where the velocity distributions are presented for a representative 
case (K = 10, m = 2, a = 0, 10, lOO) where T^<Tg. The associated 

temperature profiles are given in figure 9. Frictional heating is 
appreciable for a = 0 and 100 but not for a = 10; and although for 
a = 100 the flow is again in the same direction as for the same case 
for Tw >Ts (fig. 4(d)), the velocities are smaller for T^ < Tg. 

Unfortunately, there are no experimental results available with 
which to compare the results predicted herein. The experiments of 
Elenbaas (reference 2) were made with short plates, and, consequently, 
there were variations of velocity and temperature profiles with the 
longitudinal distance. Such variations were not considered in this 
analysis . 


Heat Transfer 

The heat-transfer coefficients for the natural-convection process 
considered here can be expressed in terms of Nusselt numbers. For the 
case where the walls are not at the same temperature (that is, m l). 



where the double subscript 0,d signifies that the temperature gradient 
is to be evaluated at either Y=0 or Ys=d, depending on which wall 
is under consideration. By means of the various transformations in the 
analysis this expression can be written 
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By use of equation (19), the Wusselt number for the wall at y = 0 is 


The Nusselt numbers were computed from these expressions over a range of 
values of E for m = -1, 0, and 2, and a = 0, 10, and 100, and are 
presented in figures 10 and 11. In general, the Nusselt numbers increase 
with increasing K, m, and a. Comparison of figures 10 and 11 shows that 
for m = -1 the Nusselt numbers for the wall at y = 1 are larger than 
for the wall at y = 0. However, as m increases to 2, this result is 
reversed . 

When the walls are at the same temperature (m = l) the Nusselt 
number can be written as 


and from the transformations in the analysis, this equation becomes 



48(2^20y + 564a + 1008m^ + 2352m + 1680) 


and for the wall at y = 1 



48 ( ' 2520 ' ) leSQm^ + 2352m + 1008) 




Again from equation (19) 



1008a 



The Nusselt numbers computed for this case are presented in figure 12. 
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CONCnJSIOKS 

An analysis of the natural-convection phenomenon shows that the 
flow and heat transfer not only are functions of the Prandtl and _ 
Grashof numhers_but also depend upon a new dimensionless factor, K. 

For values of K which are not negligible, the frictional heating and 
compression work may appreciably alter the results. Consideration was 
given to the particular simplified case of the fully developed natural- 
convection flow of fluids with and without heat sources between two long 
parallel surfaces oriented in the direction parallel to the generating 
body force. These surfaces were taken to have constant but not 
necessarily equal temperatures. The velocity and temperature distribu- 
tions for this special case were determined, and it was observed that 
increasing the wall-temperature ratio increased the flow velocities, 
the net mass through-flow, and the temperatures. The effect of the 
heat sources was also found to increase velocities, temperatures, and 
mass flows. The frictional heating appreciably altered the velocity 
and temperature profiles in some cases, showing that the velocities 
were increased and the heat transfer at the walls was greatly changed 
by this effect. Consideration of the frictional heating also led to 
the prediction of two flow and heat-transfer states for a given set of 
conditions and to a critical set of conditions beyond which no solutions 
existed'. These last two results are being more completely investigated. 

Lewis Fli^t Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, October 6, 1952 
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APPENDIX A 
SYMBOLS 

The following notation is used in this report; 
a^, n’t coefficients in successive approximation solutions 


■4 

specific heat at constant pressure 

Oy 

specific heat at constant volume 

d 

characteristic length (specifically distance "between plates) 

fi 

components of body force per unit mass, i=l, 2, 3, . . . 

fx 

negative of X-component of body force per unit mass 

Or 

PfxV^ 

Grashof number, — 

V 

g 

gravitational force per unit mass (or acceleration due to 
gravity) 

h 

heat-transfer coefficient 

K 

3fyd 

dimensionless parameter, Pr Gr 

Cp 

K 

dimensionless parameter^ x/Pr 

k 

thermal- conductivity coefficient 

m 

ratio of wall-temperature differences^ - Tg 

N 

a constant 

Nu 

Nusselt number \ 

P 

pressure 


I" - Ps 

Ps 

hydrostatic pressure 
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Pr 

P 

Q 

Ee 

T 

U 

U 

Ui 

u 

V 
^1 

Y 

y 

a 

P 

7 

0 

0 

0 * 


Prandtl number 
dimensionless pressure 
heat added by heat sources 
Eeynolds number 
temperature 

longitudinal velocity component 
characteristic reference velocity 
velocity components, i = 1, 2, 3, . . . 

prescribed reference velocity in forced -convection flow 
dimensionless longitudinal velocity component 
dimensionless velocity components, i = 1, 2, 3, . . . 

transverse velocity component 

rectangular Cartesian coordinates, i = 1, 2, 3, . . . 
dimensionless coordinates, i = 1, 2, 3, . . . 
coordinate 

dimensionless coordinate 

dimensionless heat-source parameter, Qd^/k0^ 
coefficient of volumetric expansion, p 

ratio of specific heats 

dimensionless temperature number, U^^/cp0^ 

T - Tg 

dimensionless temperature difference, 0/0^ 
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li absolute viscosity coefficient 

V kinematic viscosity coefficient 

p density 

-c dimensionless temperature 

Subscripts : 

a denotes adiabatic condition 

i, J rectangular Cartesian tensor and summation subscripts 

r successive approximation subscript 

s denotes a reference condition (usually taken as the 

hydrostatic condition) 


w 

denotes 

wall conditions 

^0 

denotes 

condition 

at 

y = 0 

Wi 

denotes 

condition 

at 

7 = 1 


Superscripts ; 

(2) denotes second flow and heat-transfer state 
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APPENDIX B 

DEEIVATION OF BUOXAECY TEEM 

It is often convenient in natural-convection studies to express 
the body force term (first term on right aide of equation (2)) aa a 
buoyancy term. To thia end, the caae ia conaidered in which the aur- 
face (apecifically here, the channel) and the fluid are at the aame 
temperature and there ia no flow. Equation (2) then becomea 

^Pa 

where the aubacript a indicates the hydrostatic condition, and 
equation (Bl) then expresses the fact that under this condition the 
body force is in equilibrium with the hydrostatic pressure gradient. 
This equilibrium is, of course, upset if there ia a temperature 
variation in the flow field, and the unbalanced force, which is the 
buoyancy force, causes a flow to be established. In order to intro- 
duce the buoyancy term into the equation, the body force and pressure 
terms in equation (2) 

Of 

^ 

can be written 

SP-t, 

Psfi + (P - Pg)fi - - -5x7 


where Pjj = P - Pg. 

Hence, in view of equation (Bl), these terms become 

SPt5 

(P - Pg)fi - (B2) 

The pressure gradient in the preceding expression appears in Just that 
form in equation (2a), but now the buoyancy term is to be expressed 
in terms of a temperature difference. It is first assumed that the 
density is a function of temperature alone, so that equation (4) C6in 
be written 


P = P(T) (B3) 

In the case of a liquid, this assumption is evident; whereas, in the 
case of a gas, it implies that the pressure changes are small as com- 
pared with the absolute pressure. If the coefficient of volumetric 
expansion P is introduced, equation (B3) can be written 
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dp = - (3pdT (B4) 

The writing of the density difference in expression (BE) in terms of 
a temperature difference can now he effected in one of several ways, 
depending on the specific problem. 

Most of the theoretical work done to date applies only for the case 
of small temperature differences. Hence, for that case the differentials 
in equation (B4) are replaced hy differences to yield 

3=1 (P~Ps) 

p (Tg-T) 

The buoyancy term in equation (BE) then becomes 

(P-Ps)fi = - PfiP(T-Tg) (B5) 

as it appears in equation (2a). 

If the analysis is not to be limited to the case of small tempera- 
ture differences, equation (B5) can be obtained in either of the 
following ways, depending on whether the fluid is a liquid or a gas. 

For liquids, it is assumed that 3P = N where N is a constant. A 
numerical check of this assumption shows that it is reasonable under 
commonly encountered conditions for most fluids; in particular, if 
the constant is evaluated at the ambient condition, the variation of 
3P over a large range of temperature is small. (For unusually large 
temperature ranges H could be evaluated at some other appropriate 
condition.) As a result of this assumption a linear density- 
temperature variation is obtained from equation (B4), and then by 
direct substitution equation (B5) is obtained. For gases, the 
equation of state is 

P = PET (B6) 

where B is the gas constant. 

Substitution into the buoyancy term yields 

^■i 

It need now only be assumed that the difference between P and P3 
is everywhere small; therefore 

Pfi 

(P-Ps)fi = (Tg-T) 


(P-Pg)fi = 


BT 


B Tg 


(B7) 
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By definition, 3 = 1/T for gases, so that the final desired form, 

(P-Pg)fi = - PfiPg(T-Tg) 

is obtained. In this case, P is evaluated at the hydrostatic condi- 
tion. An eq.uivalent form could also he obtained vhere the density 
would be evaluated at that condition and p could be variable. 
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Figure 1. - Schematic sketch of simplified configuration. 
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y = Y/d 


(a) Ratio of wall-temperature differences, m, -1. 

Figure 2, - Dimensionless velocity distributions for various heat 
source parameters and for K = 0.5. 
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Figure 2. - Continued. Dlmenelonlees velocity distributions for various 
heat-source parameters and for K = 0.5. 
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(c) Eatlo of wall-temperature differences, m, 1. 

Figure 2. - Continued. Dimensionless velocity distributions for various 
heat-source parameters and for K = 0.5. 
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Figure 2. - Concluded. Diraensionlese velocity distributions for various 
heat -source parameters and for K = 0.5. 
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(a) Eatio of wall-temperature differences^ m, -1. 

Figure 3. - Dimensionless velocity distributions for various heat-source 
parameters and for K = 3.0. 


Dimensionless velocity 


32 


NACA TN 2863 



Figure 3. - Continued. Dlmenslonlesa velocity distributions for various 
heat-source parameters and for K = 3.0. 
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(c) Ratio of wall-temperature difference, m, 1. 

Figure 3. - Continued. Dimensionless velocity distributions for various 
heat-source parameters and for K = 3.0. 
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0 .2 .4 .6 .8 1.0 

y = Y/d 


(d) Ratio of wall -temperature differences, m, 2. 

Figure 3. - Concluded. Dimensionless velocity distributions for various 
heat-source parameters and for K = 3.0, 
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(si) ■Ratio of vall-temperature differences, m, -1. 

Figure 4. - Dimensionless velocity distributions for various heat-source 
parameters and for K = 10.0. 
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(b) Katio of wall-temperature differences, m, 0. 

Figure 4, - Continued. Dimensionless velocity distributions for various 
heat-source parameters and for K = 10.0. 
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(c) Ratio of wall-temperature differences, m, 1. 

Figure 4. - Continued. Dimensionless velocity distributions for various 
heat-source parameters and for K = 10.0. 
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(d) Ratio of wall-temperature differences, m, 2, 

Figure 4. - Concluded . Dimensionless velocity distributions for various 
heat-source pareimeters for K = 10.0. 
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y = Y/d 


(a) Ratio of wall -temperature differences, m, -1. 

Figure 5, - Dimensionless temperature distribulions for various heat 
source parameters and for K = 0.5. 
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0 .2 .4 .6 .8 1.0 

y = Y/d 


(b) Ratio of wall-temperature differences, m, 0. 

Figure 5. - Continued. Dimensionless temperature distributions for various 
heat-source parameters and for K = 0.5. 
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y = Y/d 


(c) Ratio of wall-temperature differences, m, 1. 

Figure 5. - Continued. Dimensionless temperature distributions for 
various heat-source parameters and for K = 0.5. 
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(d) Ratio of wall-temperature differences, m, 2. 

Figure 5. - Concluded. Dimensionless temperature distributions for 
various heat-source parameters and for K = 0.5. 
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(a) Ratio of wall -temperature differences^ m, -1. 

Figure 6* - Dimensionless temperature distributions for various heat 
source parameters and for K = 3,0, 
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Figure 6. - Continued. Dlmensionlese temperature distributions for 
various heat-source parameters and for K = 3.0. 
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7 = Y/d 


(c) Ratio of wall-temperature differences, m, 1. 

Figure 6. - Continued. Dimensionless temperature distributions foi 
various heat-source parameters and for K = 3.0* 
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y = Y/d 

(d) Ratio of wall-temperature differences, m, 2, 

Figure 6. - Concluded. Dimensionless temperature distritutions for various 
heat -source parameters and for K = 3.0, 
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y = Y/d 

(a) Ratio of wall -temperature differences^ m, -1. 

Figure 7. - Dimensionless temperature distributions for various lieat- 
source parameters and for K = 10.0. 
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y = Y/d 

("b) Ratio of wall-temperature differences, m, 0. 

Figure 7. - Continued. Dimensionless temperature distrilDutions for various 
heat- source parameters and for K = 10.0. 
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(c) Ratio of wall-temperature differences, m, 1. 

Figure 7. - Continued. Dimensionless temperature dlstrltutlons for 
various heat -source parameters and for K = 10.0. 
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y = Y/d 

(d) Ratio of wall-temperature differences, m, 2. 


Figure 7. - Concluded. Dimensionless temperature dlstrltutlons for various 
heat -source parameters and for K = 10.0. 
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y = Y/d 

Figure 8. - Dimensionless velocity distributions for T^<Tg. K = 10.0; 
ratio of wall-temperature differences, m, 2, 
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0 .2 .4 .6 .8 1.0 

y = Y/d 


Figure 9. - Dimensionless temperature dlstrltutlons for T^<Tg. K = 10.0; 
ratio of wall -temperature differences, m, 2. 
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Figure 10. - Nusselt number for wall at y = 0 for ratio of wall-temperature 
differencee m of -1, 0, and 2. Curves for a = 100 are plotted for 
one-tenth actual value. 
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Figure 11. - Nusselt number for wall at y = 1 for ratio of wall- temperature 
differences m of -1^ 0^ and 2. Curves for a = 100 are plotted for 
one-tenth actual value. 
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Figure 12. - Nuseelt number for valla at same temperature (m = l). 


NACA-LaneUy - 12-29-52 - 1000 





